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ABSTRACT 


Methods  of  approximating  inviscid  boundaries  of  jets 
exhausting  from  axisymmetric  supersonic  nozzles  into  a 
quiescent  atmosphere  and  moving  streams  have  been 
developed.  Jet  shapes  from  two  typical  afterbody  con¬ 
figurations  were  calculated  for  various  flight  conditions. 
Good  agreement  with  available  theoretical  and  experi¬ 
mental  data  was  obtained.  A  step-by-step  calculational 
procedure  is  given  in  the  Appendix. 
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1.0  INTRODUCTION 

A  jet  exhausting  from  an  under  expanded  supersonic  nozzle  expands 
rapidly  about  the  nozzle  lip  and  forms  a  plume  of  complicated  aero¬ 
dynamic  structure.  The  boundary  of  the  plume  is  a  mixing  region 
resulting  from  the  viscous  interaction  of  the  jet  with  the  surrounding 
medivun.  In  order  to  simplify  the  problem,  the  jet  gas  and  the  external 
gas  in  the  mixing  region  may  be  assumed  to  lie  in  separate  areas  inter¬ 
nally  and  externally,  respectively,  to  a  nonviscous  streamline.  Such  an 
analogy  is  possible  since  the  total  mixing  region  acts  as  a  nonviscous 
streamline  of  finite  thickness  -  separating  the  jet  and  the  surrounding 
medium. 

A  number  of  experimental  and  theoretical  studies  have  been  reported 
on  the  spreading  of  underexpanded  jets  exhausting  into  a  quiescent  atmos¬ 
phere  based  on  the  above  analysis  of  the  jet  boundary.  However,  only  a 
few  have  considered  the  case  of  external  flow.  The  method  of  character¬ 
istics  is  a  technique  available  for  the  calculation  of  the  jet  structure. 

The  characteristics  solution  reported  by  Wang  and  Peterson  (Ref.  1)  for 
jet  spreading  in  a  hypersonic  stream  utilizes  Newton's  impact  equation 
to  determine  the  pressure  variation  along  the  jet  boundary  caused  by  the 
interaction  with  the  external  stream.  With  the  boundary  conditions 
determined,  the  method  of  characteristics  can  be  readily  applied  to 
calculate  the  jet  flow  field.  However,  the  computations  are  lengthy  and 
require  electronic  computing  equipment.  To  avoid  the  complications  of 
the  characteristics  solution,  Latvala  (Ref.  2)  and  Adamson  and  Nicholls 
(Ref.  3)  replaced  the  effects  of  expansion  waves  impinging  on  the  jet 
boxmdary  by  a  quasi  one -dimensional  area  increase,  and  developed  a 
simplified  method  of  estimating  jet  spreading  into  a  quiescent  atmos¬ 
phere. 

The  present  study  adapts  Newton's  impact  theory  as  in  Wang  and 
Peterson's  solution  and  the  conical  flow  theory  as  applied  to  the  flow 
about  bodies  of  revolution  to  a  one -dimensional  flow  solution.  The  results 
are  simple  approximate  solutions  for  estimating  the  boundaries  of  jets 
expanding  into  either  a  quiescent  atmosphere  or  an  external  flow. 

Satisfactory  comparison  with  theoretical  results  obtained  by  the 
characteristics  solution  and  Latvala 's  solution  are  shown  for  the  quiescent 
atmosphere  condition.  Only  a  limited  amount  of  jet  spreading  data  with 
external  flow  was  available  for  comparison  with  the  methods  suggested 
here;  however,  satisfactory  comparisons  were  made  which  should  justify 
the  use  of  the  methods  for  further  studies. 


Manuscript  released  for  publication  March  1962. 
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2.0  ANALYSIS 


The  present  analysis  of  a  supersonic  jet  expanding  into  a  quiescent 
atnaosphere  or  a  high  velocity  stream  employs  some  of  the  principal 
assumptions  and  equations  used  by  Wang  and  Peterson  (Ref.  1)  and 
Latvala  (Ref.  2) 

In  the  case  of  an  underexpanded  nozzle,  there  are  expansion  waves 
originating  at  the  nozzle  lip  and  propagating  downstream  toward  the  jet 
boundary.  These  waves  impinge  on  the  boundary,  turn  the  flow,  and  are 
reflected  as  compression  waves.  The  compression  waves  then  converge 
to  form  a  finite  shock  wave.  The  method  of  characteristics  reported  by 
Wang  and  Peterson  is  based  on  this  analysis  of  the  jet  structure. 

The  present  methods  and  the  method  developed  by  Latvala  depend 
on  replacing  the  effects  of  expansion  waves  impinging  on  the  jet  boimdary 
flow  by  one -dimensional  flow  relations.  However,  the  present  method 
is  different  from  Latvala 's  method  in  that  Newtonian  theory  is  used  in 
conjunction  with  one -dimensional  flow  theory  to  define  the  jet  structure. 

By  one -dimensional  theory,  as  the  area  of  the  expanding  jet  increases 
downstream  of  the  nozzle  exit,  the  internal  boundary  pressure  decreases. 
To  maintain  boundary  equilibrium,  the  boundary  flow  direction  changes, 
retarding  the  expansion  of  the  jet.  Compression  waves  generated  from 
the  boundary  converge  to  form  a  shock  originating  immediately  down¬ 
stream  of  the  nozzle  lip.  The  shock  structure  is  analogous  to  the  struc¬ 
ture  described  above  for  the  method  of  characteristics.  For  jet  boimdary 
Mach  numbers  in  the  hypersonic  range,  the  boundary  shock  is  close  to 
the  boundary  streamline  (see  photographs  in  Ref.  2).  Therefore,  there 
will  be  very  little  flow  across  the  shock,  and  the  Newtonian  impact  theory 
is  justified  for  the  internal  jet  flow  as  is  the  assumption  of  an  isentropic 
flow  field.  In  the  present  analysis,  only  the  case  where  the  internal  flow 
near  the  boundary  was  in  or  near  the  hypersonic  range  (Mb  ^  4.  0)  was 
considered.  However,  for  the  case  where  the  flow  near  the  boundary  is 
in  the  supersonic  range,  it  may  be  possible  to  incorporate  oblique  shock 
relations. 

Two  typical  afterbody  configurations  are  considered.  For  Configura¬ 
tion  1  the  nozzle  exit  diameter  is  equal  to  the  model  base  diameter 
(Fig.  la),  whereas  for  Configuration  2  the  nozzle  exit  diameter  is  smaller 
than  the  model  base  diameter  (Fig.  lb).  Although  there  are  many  different 
base  configurations,  the  most  typical  may  be  grouped  into  these  two  basic 
categories.  The  afterbody  configuration  does  not  influence  the  boundary 
of  a  jet  expanding  into  a  quiescent  atmosphere;  however,  its  influence  is 
appreciable  for  external  flow. 
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2.1  QUIESCENT  ATMOSPHERIC  CONDITIONS 

When  a  jet  is  exhausting  into  a  quiescent  atmosphere,  it  may  be 
assumed  that  the  pressure  along  the  external  boundary  remains  constant 
and  is  equal  to  that  of  the  ambient  conditions.  By  one -dimensional  flow 
relations,  the  internal  flow  will  vary  between  points  on  the  boundary. 
Consequently,  only  one  set  of  flow  equations  representing  the  internal  jet 
conditions  is  necessary  to  describe  the  jet  boundary. 

The  initial  expansion  of  the  jet  will  follow  Prandtl-Meyer  flow  as 
outlined  in  Refs.  2  and  3.  The  expansion  is  complete  when  the  pressure 
on  the  jet  boundary  has  reached  the  ambient  conditions.  Prandtl-Meyer 
equations  may  be  used  to  represent  the  initial  expansion  of  the  internal 
flow 

Oi  =  I'l  -  +  0N  (1) 

where 


The  notation  employed  in  Eq.  (1)  is  indicated  in  Fig.  la. 

The  shape  of  the  jet  boundary  downstream  of  the  nozzle  exit  may  be 
approximated  by  considering  quasi  one -dimensional  area  increases  of 
the  jet.  That  is,  it  may  be  assumed  that  one  can  find  the  Mach  number 
and  pressure  at  any  axial  position  in  an  expanding  flow  from  the  area 
ratio  at  a  given  point.  Therefore,  (A/A*)i  corresponding  to  (p^/Pq)  can 
be  obtained  from  compressible  flow  tables  based  on  the  specific  heat 
ratio  of  the  exhausting  jet.  By  beginning  with  a\  and  (A/A*)i  at  the  noz¬ 
zle  lip,  a  point -by-point  calculation  along  the  jet  boundary  can  be  made. 

To  determine  the  jet  boundary  shape  it  is  necessary  to  assume  a 
one -dimensional  area  increase  of  the  jet  over  a  finite  distance  between  a 
point  on  the  boundary  denoted  by  (n)  at  a  radius  (r^j)  and  a  nearby  point 
(n  +  1)  at  a  radius  rn  +  Ar.  Flow  conditions  within  the  jet  are  assumed 
to  be  constant  over  a  spherical  surface,  illustrated  on  the  following  page 
as  a  circular  arc,  centered  at  the  hypothetical  source  A. 
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The  smaller  the  value  used  for  Ar  the  better  will  be  the  approximation. 
With  Ofn  calculated  and  Ar  assumed  between  (n)  and  (n  +  1)  the  boundary- 
point  (n  +  1)  may  be  located  by 

+  I  — 

'‘i  "  "j  ' 

+  I  „ 

The  external  ambient  pressure  is  constant  between  (n)  and  (n  +  1), 
whereas  the  internal  pressure  has  decreased  corresponding  to  one¬ 
dimensional  flow  relations.  At  (n  +  1)  the  boundary  is  no  longer  in  equi¬ 
librium.  To  balance  the  pressure  internally  and  externally,  it  is  neces¬ 
sary  to  decrease  the  angle  of  the  boundary  at  (n  +  1),  compressing  the 
flow  near  the  boundary  to  the  ambient  pressure.  With  the  resulting 
angle  Qfn+  1  rns^y  define  a  new  hypothetical  source  B  as  shown  in  the 
sketch.  The  jet  flow  field  is  then  assumed  to  be  related  to  the  ambient 
pressure  which  exists  both  internally  and  externally  to  the  jet  boundary 
immediately  aft  of  the  point  (n  +  1).  Beginning  with  the  new  flow  condi¬ 
tions  a'n+  1  and  (A/A*)  (corresponding  to  Pa,/Pc)  one  may  now  assume  a 
one -dimensional  area  increase  from  (n  +  1)  to  (n  +  2)  with  B  as  the  hypo¬ 
thetical  source.  This  analysis  is  similar  to  that  used  at  the  nozzle  exit. 
As  the  jet  expands  about  the  nozzle  lip  to  the  ambient  pressure,  the 
value  of  (A /A*)  increases  correspondingly.  During  the  expansion  the 
hypothetical  source  moves  from  P  to  Q. 

To  determine  the  resulting  angle  an+  1  the  assumption  is  made  that 
the  varying  internal  flow  between  (n)  and  (n  +  1)  may  be  represented  by 


‘j 

Ax 


Ar  cot  On 


(2) 

(3) 
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average  flow  conditions  between  these  two  points  if  the  corresponding 
radius  change  is  small.  This  approximation  is  analogous  to  that 
employed  in  the  generalized  shock  expansion  method.  Ref.  4,  where 
curved  body  shapes  are  approximated  by  a  series  of  conical  frustums. 

In  this  analysis  it  allows  the  flow  at  the  point  (n)  to  be  simply  related  to 
the  flow  at  point  (n  +  1).  Since  the  jet  flow  of  a  highly  underexpanded 
nozzle  is  in  or  near  the  hypersonic  range  (Mb  2  4.  0),  the  Newtonian 
impact  equation  may  be  applied  to  the  internal  flow.  Employing  New¬ 
tonian  theory  and  representing  the  flow  near  the  boundary  over  a  small 
interval  by  the  average  conditions  in  the  interval  results  in  the  following 
equation. 


For  the  internal  flow  the  dimensionless  Newtonian  impact  equation 


IS 


Pb 


n  +  1 


(Pb). 


(“n  -  “n  +  l)  +  l] 


(4) 


where  Pbn+  i  “  P®  quiescent  atmosphere  case. 


From  Eq.  (4)  the  new  angle  is 


Po. 


(Pb) 


avg 


(pb) 


avg 


7:  (Mb) 

J  avg 


(4a) 


The  average  internal  flow  conditions  between  (n)  and  (n  +  1)  corre¬ 
spond  to  the  average  Mach  number  between  (n)  and  (n  +  1),  assuming 
the  Mach  number  varies  linearly  with  jet  radius.  Mach  number  was 
used  as  a  basis  for  the  average  flow  conditions  rather  than  pressure 
since  the  change  in  Mach  number  for  change  in  jet  radius  is  more  nearly 
linear  (see  Fig.  2).  This  assumption  is  reasonable  for  small  Ar.  The 
average  Mach  number  may  be  calculated  from 


(Mb) 


avg 


Mbp  +  ^bp  -)-  I 
2 


(5) 


Mbn  and  Mbn+  i  correspond  to  (A/A*)j^  and  (A/A*)jj+  (A/A*)^  is 

known  and  {A/A*)^^  may  be  calculated  as  follows: 


{rny 


(6) 


From  compressible  flow  tables  based  on  yj,  (Pb)avg^Pc  correspond¬ 
ing  to  (Mb)^yg  may  be  obtained. 
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An  approximate  nearby  boundary  point  (n  +  2)  may  now  be  located 
by  assuming  the  boundary  flow  angle,  0?^+  1»  constant  from  (n  +  1) 

to  (n  +  2).  The  location  is  given  by  Eqs.  (2)  and  (3).  This  procedure  is 
begun  at  the  nozzle  lip  (n  =  1)  and  repeated  application  will  provide  a 
step-by-step  estimation  of  the  jet  boundary. 

To  facilitate  the  use  of  this  method,  a  step-by-step  procedure  is 
given  in  the  Appendix. 

A  mathematically  correct  solution  for  calculating  the  jet  boundary 
would  be  to  apply  the  correction  angle  required  to  place  the  inviscid 
streamline  in  pressure  equilibrium  at  (n  +  1)  to  the  previous  point  (n). 
However,  it  was  found  that  this  solution  produced  a  more  compressed 
jet  boundary  than  actually  exists.  It  is  shown  later,  in  comparisons 
with  experimental  boundaries,  that  the  presented  methods  result  in  very 
good  agreement.  An  explanation  for  the  greater  expansion  of  the  experi¬ 
mental  jet  may  be  credited  to  the  initial  overexpansion  of  the  jet  about 
the  nozzle  lip.  Also,  for  the  case  of  external  flow,  the  compression 
caused  by  the  free  stream  impinging  on  the  jet  was  assumed  to  be  instan¬ 
taneous  at  a  point;  whereas,  the  compression  actually  takes  place  over 
a  small  interval. 

It  is  interesting  to  note  that  the  final  results  of  applying  the  cor¬ 
rected  streamline  angle  at  (n  +  1)  instead  of  (n)  will  differ  from  zero  at 
the  nozzle  lip  to  the  Ar  value  at  the  maximum  jet  diameter  (or  =  0).  This 
indicates  that  the  two  methods  are  essentially  equivalent  for  sufficiently 
small  values  of  Ar. 


2.2  HYPERSONIC  EXTERNAL  FLOW -CONFIGURATION  1 

The  shape  of  a  jet  boundary  which  is  exposed  to  an  external  stream 
is  influenced  by  the  flow  conditions  both  internal  and  external  to  the  jet. 
Consequently,  two  sets  of  flow  equations,  one  representing  the  internal 
jet  conditions  and  the  other  representing  the  external  flow,  are  necessary 
to  describe  the  jet  boundary.  Variables  satisfying  both  the  internal  and 
external  flow  equations  simultaneously  specify  the  boundary. 

Similar  to  the  expansion  of  a  jet  into  a  quiescent  atmosphere,  the 
initial  expansion  of  the  jet  into  an  external  stream  will  follow  Prandtl- 
Meyer  flow.  The  expansion  is  complete  when  the  pressure  on  the  jet 
boundary  is  in  equilibrixun  with  the  surrounding  medium.  As  before, 

Eq.  (1)  may  be  used  to  represent  the  initial  expansion  of  the  internal 
flow: 

=  *^1  -  V.  + 
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When  a  jet  is  issuing  into  a  moving  stream,  the  boundary  pressure 
is  no  longer  constant  and  the  interaction  between  the  jet  and  the  external 
stream  must  be  taken  into  accoimt.  As  far  as  the  external  flow  is  con¬ 
cerned,  the  jet  may  be  treated  as  a  solid  body,  and  with  a  hypersonic 
external  flow,  Newtonian  theory  can  be  used.  The  dimensionless  New¬ 
tonian  impact  equation  representing  the  external  boundary  conditions  at 
the  nozzle  lip  (n  =  1)  is  defined  by 

[y*  sin*  +  1]  (7) 

The  initial  expansion  equation  and  Eq.  (7)  are  both  functions  of  aj 
and  Pb]^/Pc’  Upon  leaving  the  nozzle  exit,  the  flow  about  the  lip  will 
expand  to  the  angle  ai  and  to  the  boundary  pressure,  p^j^.  It  is  evident 
that  only  one  pair  of  values  of  ai  and  Pbj^/Pc  define  the  initial  portion 
of  the  boundary.  These  values  may  be  determined  by  simultaneously 
solving  Eqs.  (1)  and  (7),  either  analytically  or  graphically,  as  shown  in 
Fig.  3  for  the  flow  conditions  used  by  Wang  and  Peterson  (Ref.  1). 


The  boundary  shape  is  determined  by  an  analysis  similar  to  that 
used  for  the  quiescent  atmospheric  condition.  The  boundary  flow  direc¬ 
tion  is  assumed  to  remain  constant  from  (n)  to  (n  +  1)  and,  knowing  orn 
and  assuming  some  Ar,  the  nearby  boundary  point  (n  +  1)  is  defined  by 
Eqs.  (2)  and  (3).  To  locate  the  boundary  point  (n  +  2),  the  boundary  angle 
at  (n  +  1)  must  be  calculated.  This  is  done  by  again  assuming  that  the 
varying  internal  flow  between  (n)  and  (n  +  1)  may  be  represented  by  the 
average  flow  conditions  between  these  two  points.  The  average  conditions 
corresponding  to  the  average  Mach  number  between  (n)  and  (n  +  1)  are 
calculated  as  described  for  the  quiescent  external  conditions.  Since  the 
external  Mach  number,  M„,  and  the  jet  boundary  Mach  number,  (Mb)a„o.< 
are  both  assumed  in  the  hypersonic  range,  the  dimensionless  Newtonian 
impact  equation  (Eq.  (7) )  may  be  used  to  determine  the  flow  conditions 
necessary  to  restore  the  external  and  internal  flow  to  equilibrium  at 
(n  +  1). 


For  the  internal  flow,  Eq.  (4)  is  used: 


>n  +  1 


(Pb) 


avg 


[yj  +  +  1] 
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whereas  for  the  external  flow,  Eq.  (7)  is  used: 


P'’  n  +1 

Pc 


sin’  «n  +  1  +  l] 


Free-stream  pressure  is  determined  from  external  flow  conditions. 

From  the  two  equations  the  angle  of  flow,  an+  1*  obtained  by 

graphical  methods.  By  again  assuming  one-dimensional  flow  between 
(n  +  1)  and  (n  +  2)  for  some  Ar,  the  boundary  point  (n  +  2)  may  be  located 
by  Eqs.  (2)  and  (3).  This  procedure  is  begun  at  the  nozzle  lip  and  is 
repeated  at  each  boundary  point. 


2.3  HYPERSONIC  EXTERNAL  FLOW -CONFIGURATION  2 

In  analyzing  the  jet  spreading  for  Configuration  2  (see  Fig.  lb),  the 
free  stream  was  assumed  to  remain  parallel  to  the  jet  axis  from  the 
missile  base  to  the  point  of  intersection  with  the  jet  boundary.  With  this 
assumption,  the  base  pressure  would  be  equal  to  the  free-stream  pres¬ 
sure  and  remain  constant  along  the  boimdary  from  the  nozzle  lip  to  the 
point  of  intersection.  If,  however,  the  missile  base  pressure  is  known 
or  can  be  determined,  the  turning  of  the  free  stream  at  the  base  may  be 
calculated  and  the  point  of  intersection  accurately  located.  At  hyper¬ 
sonic  velocities  the  base  pressure  approaches  the  free-stream  pressure, 
and  subsequent  turning  of  the  free  stream  would  be  small. 

The  initial  portion  of  the  boundary  from  the  nozzle  lip  to  the  point 
of  intersection  with  the  external  stream  is  assumed  to  be  a  boundary 
with  constant  external  pressure  and  may  be  determined  as  described  for 
the  quiescent  condition.  At  the  point  of  impact  the  external  stream  begins 
to  retard  the  jet  expansion  and  the  boundary  pressure  is  no  longer  con¬ 
stant.  It  may  then  be  assumed  that  the  flow  conditions  at  the  point  of 
intersection  are  similar  to  the  flow  conditions  about  the  lip  of  an  over- 
expanded  hypersonic  nozzle,  as  illustrated  below. 
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If  the  external  Mach  number,  M„,  and  the  jet  Mach  number,  (Mb)avg» 
at  the  point  of  intersection  are  both  in  the  hypersonic  range,  the  impact 
equation  may  be  used  to  determine  the  boundary  pressure  caused  by 
internal  and  external  flow  at  the  intersection  of  both  streams. 


For  the  internal  flow,  Eq.  (4)  is  applicable  (n  =  i-1) 


•  *  ! 

Sin  (  a .  _ 


Ml, .  +  Mu .  (  Pk  ) 

where  (Mb)  =  — 1 — — 5—  and  - corresponds  to  (Mb) 

8Vg  2  Pj,  avg 

(See  the  analysis  on  average  internal  flow  conditions  as  described  in 
section  2.1.) 


For  the  external  flow,  Eq.  (7)  is  used: 


Pbi 


Pr 


fy  M  *  sin*  «•  +  ll 

j_  '  oo  00  1  J 


These  two  equations  specify  the  angle  of  the  flow,  ai,  and  the  pressure, 
P|^,,  at  the  intersection  point  of  the  jet  and  the  free  stream.  The  new 
value  of  M^j  can  then  be  obtained  as  a  function  of  p^^^/Pc  Calculations 
of  subsequent  points  are  the  same  as  for  Configuration  1  using  the  condi¬ 
tions  calculated  after  intersection  has  occurred. 

To  facilitate  the  use  of  this  method  for  both  configurations,  a  step- 
by-step  procedure  is  given  in  the  Appendix. 


2.4  SUPERSONIC  EXTERNAL  FLOW 

When  a  jet  is  issuing  into  a  supersonic  stream,  the  external  bound¬ 
ary  pressure  may  no  longer  be  defined  by  Newtonian  theory.  The  present 
approximation  involves  compressing  the  free  stream  through  a  conical 
shock  created  by  the  interaction  of  the  free  stream  with  the  jet  at  the 
nozzle  lip.  Aft  of  the  shock,  the  external  flow  is  defined  by  a  series  of 
small  angle  Prandtl-Meyer  expansions  about  the  jet  boundary. 


The  pressure  coefficient  equation  is 

(8) 

where  Cp  =  f(Qri,  M„)  (defined  from  conical  flow  charts). 


Cp  = 


Pb,-P„ 


y  M  ’ 

w  ou 


-  1 
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Rearranging  Eq.  (8)  the  initial  boundary  pressure  ratio  may  be 
expressed  as  a  function  of  the  free-stream  conditions: 


y  M  ’ 

oo  .00  p  , 

2  ^ 


(8a) 


The  initial  expansion  of  the  jet  follows  Prandt -Meyer  flow,  as  pre¬ 
viously  discussed.  Unique  values  of  the  initial  flow  conditions  (ai  and 
Pbi/Pc)  may  be  determined  by  simultaneously  solving  Eqs.  (1)  and  (8a) 
either  analytically  or  graphically. 


The  boundary  shape  is  determined  by  an  analysis  similar  to  that 
used  for  the  hypersonic  external  flow  case.  The  difference  is  in  the 
definition  of  the  external  boundary  pressure. 


For  the  external  flow, 

X  1  “  ^  +  1) 

n  +  1  n  '  ' 


(9) 


where 


Pb 


n  +  1 


=  f  (I'. 


and 


Pb 


n  +  1 


IW.  (^)(^) 

Pt  We/VPb„/ 


For  the  internal  flow,  Eq.  (4)  is  used: 

[>'i  (“i  -  “n  +  l)  +  ‘] 


From  the  above  equations,  an  incremental  analysis  of  the  boundary 
is  made,  beginning  with  the  initial  conditions  at  the  nozzle  lip. 

The  analysis  of  the  jet  boundary  for  Configuration  2  is  also  similar 
to  the  corresponding  analysis  for  hypersonic  external  flow.  The  quies¬ 
cent  calculations  are  used  to  determine  the  initial  portion  of  the  botmd- 
ary  to  the  point  of  intersection  with  the  external  stream.  At  the  inter¬ 
section,  assuming  an  overexpanded  nozzle  to  exist,  the  supersonic  flow 
method  as  described  in  this  section  is  incorporated  to  continue  along  the 
boundary. 

To  facilitate  the  use  of  this  method  for  both  configurations,  a  step- 
by-step  procedure  is  given  in  the  Appendix. 


3.0  EXPERIMENTAL  INVESTIGATION 

Experimental  data  used  in  the  present  study  were  obtained  from 
wind  tunnel  tests  conducted  by  the  Research  and  Performance  Branch, 


10 


AEDC.TDR-62-7 


Propulsion  Wind  Tunnel  Facility  (PWT),  in  the  12-in.  Supersonic  Tun¬ 
nel  of  the  von  Karman  Gas  Dynamics  Facility  (VKF)  at  the  Arnold 
Engineering  Development  Center  (AEDC). 


3.1  TUNNEL  DESCRIPTION 

The  12-in.  Supersonic  Tunnel  is  an  intermittent  wind  tunnel  with  a 
Mach  number  range  from  1.5  to  5.  Stagnation  pressures  are  regulated 
with  automatic  control  valves  which  throttle  the  flow  from  the  air- 
storage  system.  Low  discharge  pressures  are  maintained  by  means  of 
a  vacuum  sphere.  For  a  more  detailed  description  of  the  12-in.  Super¬ 
sonic  Tunnel  and  its  operating  parameters  see  Ref.  5. 


3.2  MODEL  DESCRIPTION 

Two  high  pressure,  cold-flow  model  configurations  were  investi¬ 
gated  (see  Fig.  4).  A  model  of  boattail  design  representing  a  typical 
no-base  afterbody  is  shown  in  Fig.  4a.  The  model  consisted  of  a 
2.  0-in.  -diam  stainless-steel  cylinder  to  which  an  ogive  nose  and  inter¬ 
changeable  exit  nozzles  could  be  attached.  In  the  present  program  only 
the  complete,  boattailed,  2.  5  exit  Mach  number  nozzle  was  investigated. 

The  model  simulating  an  afterbody  configuration  with  a  model  base 
diameter  larger  than  the  nozzle  exit  diameter  (see  Fig.  4b)  consisted  of 
a  2.  0-in.  stainless -steel  cone-cylinder  to  which  interchangeable  nozzles 
could  be  attached.  In  the  present  study  only  the  3.  7  exit  Mach  number 
nozzle  was  investigated. 


3.3  DATA  ANALYSIS 

Experimental  data  used  in  comparison  with  the  present  method  were 
obtained  by  the  schlieren  photographic  and  shadowgraph  processes.  The 
knife  edge  for  schlieren  photographs  in  this  report  was  oriented  such  that 
positive  density  gradients  in  a  vertical  downward  direction  appear  as  dark 
regions  on  the  photograph. 

A  typical  schlieren  photograph  of  a  jet  issuing  from  a  supersonic 
nozzle  into  a  hypersonic  free  stream  is  given  in  Fig.  5.  The  boundary 
of  the  jet  is  defined  by  a  mixing  region,  as  noted  on  the  photograph.  The 

doncity  gradion-t  of  tho  ziaixing  boundary  normal  to  tbo  knifo  edge  is  posi¬ 
tive,  as  indicated  by  the  darkness  of  the  photograph. 
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The  free  stream  impinging  on  the  mixing  boundary  is  initially  com¬ 
pressed  through  a  conical  shock  and  then  expands  about  the  boundary 
similar  to  flow  about  a  body  of  revolution.  The  density  gradient  across 
the  shock  normal  to  the  knife  edge  is  positive  and  appears  as  a  dark 
region  on  the  schlieren  photograph.  The  expansion  of  the  free  stream 
aft  of  the  shock  creates  a  negative  density  gradient  normal  to  the  knife 
edge  and  appears  as  a  light  region  on  the  schlieren  photograph. 

Internal  to  the  mixing  boundary  is  a  shock  formed  by  expansion 
waves  from  the  nozzle  reflecting  from  the  boundary  as  compression 
waves.  Jet  flow  passing  through  the  shock  is  compressed,  creating  a 
negative  density  gradient  normal  to  the  knife  edge.  The  corresponding 
area  in  Fig.  5  is  light,  as  noted. 

A  shadowgraph  of  the  jet  structure  may  be  analyzed  similar  to  that 
above,  noting  that  a  shadowgraph  measures  rate  of  change  of  density 
gradient. 


4.0  DISCUSSION 


4.1  COMPARISON  WITH  THEORY 

Boundaries  calculated  by  the  present  method  for  quiescent  atmos¬ 
pheric  conditions  are  compared  with  theoretical  boundaries  obtained  by 
the  method  of  characteristics  and  Latvala's  solution  (results  taken  from 
Ref.  2)  in  Fig.  6.  Close  agreement  was  obtained  at  pressure  ratios 
(Pj/Poo)  of  81  and  above.  A  deviation  of  boundaries  was  obtained  at  a 

pressure  ratio  (pj/Po,)  of  10,  with  the  present  solution  giving  the  smaller 
boundary. 

Figure  7  presents  a  comparison  of  the  present  hypersonic  method 
with  the  characteristics  solution  by  Wang  and  Peterson  at  free-stream 
Mach  numbers  of  4,  5,  and  6.  Boundaries  calculated  by  the  hypersonic 
method  give  about  15  percent  less  expansion  in  radius  than  the  character¬ 
istics  solution.  The  pressure  at  points  on  the  boundary  caused  by 
external  flow  is  calculated  using  Newton's  impact  equation  in  both  methods; 
however,  computations  to  correct  the  internal  conditions  for  the  change 
in  external  pressure  are  considerably  different,  as  described  in  the 
analysis. 

4.2  COMPARISON  WITH  EXPERIMENTAL  DATA 

Jet  boundaries  calculated  by  the  present  method  and  in  Refs.  1 
thr.ough  3  are  for  inviscid  flow  and  essentially  represent  dividing  stream¬ 
lines  between  the  jet  and  external  stream.  When  these  inviscid  boundaries 
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are  compared  with  experimental  boimdaries,  they  should  obviously  lie 
within  the  viscous  mixing  region  (see  Fig.  5);  however,  their  exact 
relation  to  the  mixing  region  is  not  known  unless  some  characteristic 
mixing  profile  is  measured.  In  order  to  relate  computed  boundaries 
to  physical  viscous  flows,  the  nonviscous  streamlines  have  been  used  to 
compute  theoretical  jet  shocks  which  then  can  be  directly  compared  with 
experimental  jet  shocks.  The  method  of  estimating  the  shock  location 
is  contained  in  Ref.  6.  It  consists  of  determining  the  local  Mach  line  at 
each  calculated  boundary  point  for  the  jet  flow.  When  these  Mach  lines 
are  graphically  located,  they  coalesce  into  the  theoretical  jet  shock. 

Boundaries  with  which  to  compare  the  present  hypersonic  method 
for  Configuration  1  were  obtained  with  a  model  of  boattail  design  (see 
Fig.  4a).  The  external  flow  was  assumed  to  be  equal  to  free-stream 
tunnel  conditions  on  the  cylindrical  portion  of  the  model,  and  then  ex¬ 
panded  (by  Prandtl-Meyer  theory)  through  the  angle  of  the  boattail  prior 
to  impinging  on  the  jet.  Calculations  were  then  made  as  described  for 
Configuration  1.  Experimental  data  for  external  Mach  numbers  of  4  and 
5  were  taken  from  shadowgraphs.  Because  of  the  poor  quality  of  the 
shadowgraphs,  only  the  jet  shocks  could  be  identified.  Boundaries  were 
calciilated  by  the  hypersonic  method,  and  from  the  boundaries  theoretical 
jet  shocks  were  constructed.  The  present  solutions  are  compared  with 
the  experimental  jet  shocks  in  Fig.  8.  Excellent  agreement  with  the  jet 
shocks  in  Fig.  8  seem  to  indicate  that  the  calculated  nonviscous  stream¬ 
lines  are  correctly  located  in  the  boundary  mixing  region.  The  initial 
portion  of  the  jet  shocks  in  Fig.  8a  could  not  be  clearly  identified  for  an 
axial  distance  of  approximately  two  nozzle  diameters  aft  of  the  nozzle 
exit.  However,  good  agreement  over  the  identifiable  portion  was  ob¬ 
tained. 

Experimental  data  to  compare  both  the  hypersonic  and  supersonic 
methods  for  Configuration  2  were  obtained  from  schlieren  photographs. 
Base  pressures  from  four  orifices  located  in  the  model  base  region 
were  used  to  correctly  locate  the  interaction  of  the  jet  boundary  and  the 
free  stream.  The  pressure  readings  were  averaged  and  the  external 
flow  -  assumed  to  be  equal  to  free-stream  tunnel  conditions  on  the  model 
afterbody  -  was  expanded  or  compressed  to  satisfy  base  conditions. 
Calculations  to  determine  the  boundary  streamline  were  then  made  as 
described  in  the  analysis.  Boundaries  calculated  by  the  present  methods 
and  the  corresponding  jet  shocks  are  compared  with  experimental  jet 
structures  in  Figs.  9,  10,  and  11  for  various  flight  conditions.  Good 
agreement  with  the  shock  was  obtained  over  the  initial  portion  of  the  jet. 

li  should  bo  notod  ibai  a.fi  of  iho  2x1.c1xixn.uzx1  diazxioicz’  of  fbc  c»poz*izxioziial 

jet  shock,  the  calculated  shock  does  not  agree  as  well  with  experiment 
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(see  Figs.  9  and  10).  This  indicates  that  the  assumption  of  one¬ 
dimensional  flow  field  becomes  less  valid  with  increasing  divergence 
of  the  jet  shock  and  the  jet  boundary. 

Comparisons  of  both  the  hypersonic  and  supersonic  methods  with 
experimental  boundaries  were  made  at  free -stream  Mach  numbers  of 
3.5,  4.0,  and  5.  0  (see  Figs.  9b,  10b,  and  11b).  Very  close  agreement 
of  the  two  methods  was  obtained  for  this  Mach  number  range.  Additional 
comparisons  of  boundaries  for  free-stream  Mach  numbers  of  2.  5  and 
8.  0  were  made  (comparisons  not  shown  in  report).  The  result  was  a 
deviation  of  7.  4  percent  and  6.  7  percent  of  the  calculated  boundary  radii 
at  4  jet  exit  radii  downstream  of  the  jet  exit  for  the  two  Mach  numbers, 

2.  5  and  8.  0,  respectively.  For  both  cases  the  hypersonic  method  pre¬ 
dicted  the  larger  boundary. 


4.3  EFFECTS  OF  VARIOUS  PARAMETERS 

The  present  methods  were  also  used  to  study  the  effects  of  various 
parameters  on  jet  spreading.  Comparisons  of  calculated  boundaries  for 
the  two  typical  afterbody  configurations  are  presented  in  Fig.  12.  As 
the  base -to -nozzle  exit  diameter  ratio  increases,  the  jet  spreading  also 
increases.  The  effects  of  external  Mach  number,  nozzle  Mach  number, 
and  static  pressure  ratio,  pj/p*,  can  be  seen  in  Figs.  13,  14,  and  15. 

It  is  apparent  that  jet  spreading  increases  with  either  an  increase  in 
Pj/Pco  or  a  decrease  in  nozzle  Mach  number,  and  decreases  with  an  in¬ 
crease  in  external  Mach  number. 

At  constant  values  of  nozzle  area  ratio,  chamber  total  pressure, 
and  altitude,  as  nozzle  Mach  number  varies  with  specific  heat  ratio, 
the  jet  boundary  will  follow  a  definite  pattern  of  change  (see  Figs.  16 
and  17).  The  relative  change  in  spreading  is  indicated  by  examination 
of  the  plot  of  the  Prandtl-Meyer  and  Newtonian  equations  in  Fig.  16.  As 
the  specific  heat  ratio  decreases,  the  initial  angle  increases  for  a  con¬ 
stant  free-stream  Mach  number.  The  jet  spreading  also  increases,  as 
shown  for  an  external  Mach  number  of  5.  0  in  Fig.  17. 


4.4  LIMITATIONS 

It  shoiald  be  pointed  out  that  the  present  methods  have  certain  in¬ 
herent  limitations  because  of  the  assumption  of  a  one-dimensional  jet 

flow  field.  Ill  CLdditioxi  io  iKio  iniiiol  Q.oourxxpiion.,  iho  folio-wing  lizniio.- 

tions  exist:  First,  the  present  solutions  do  not  give  a  means  of  calcu¬ 
lating  the  mixing  region  of  the  jet  boundary.  Second,  the  supersonic 
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external  flow  method  cannot  be  used  when  the  external  flow  separates 
from  the  model  afterbody  prior  to  impinging  on  the  jet  boundary,  unless 
the  point  and  angle  of  separation  are  known  or  may  be  determined.  Third, 
the  boundary  Mach  number  must  be  hypersonic  in  order  to  use  Newtonian 
theory  for  the  internal  flow  of  the  jet.  Fourth,  the  solutions  are  based  on 
a  continuum  theory;  and  lastly,  the  methods  are  only  valid  when  the  un¬ 
disturbed  free -stream  flow  is  axisymmetric  to  the  jet  axis  (model  angle 
of  attack  of  zero).  It  should  be  noted  that  with  the  exception  of  the  first 
and  fourth  limitations,  all  concern  special  cases.  Under  practical 
operating  conditions,  the  methods  presented  will  produce  a  good  approxi¬ 
mation  of  the  jet  spreading. 

Although  the  present  method  deals  with  equilibrium  flow  within  the 
jet,  it  may  be  noted  that  a  modification  of  the  present  method  may  pos¬ 
sibly  be  used  when  the  jet  flow  is  not  in  equilibrium.  If  the  7j  distribu¬ 
tion  along  the  axis  of  the  jet  is  known,  or  may  be  determined,  the  calcu¬ 
lation  over  each  increment  on  the  boundary  may  be  based  on  the 
corresponding  7  j.  This  technique  would  be  based  on  equilibrium  flow 
only  over  small  increments  of  the  jet. 


5.0  CONCLUDING  REMARKS 


The  present  method,  as  used  for  quiescent  external  conditions, 
shows  good  agreement  with  the  method  of  characteristics  and  the  method 
introduced  by  Latvala  (Ref.  2).  A  comparison  of  the  present  method 
with  experimental  data  for  a  quiescent  atmosphere  was  not  deemed 
necessary  since  both  the  method  by  Latvala  and  the  characteristics 
solution  have  previously  shown  good  agreement  with  experimental  data 
for  the  quiescent  external  condition. 

A  comparison  of  the  supersonic  and  hypersonic  methods  with  the 
characteristics  solution  for  external  Mach  numbers  of  4,  5,  and  6  shows 
slightly  lower  boundaries,  the  difference  in  boundary  diameter  being  less 
significant  at  the  larger  free-stream  Mach  numbers.  However,  the 
present  method  results  for  the  case  of  supersonic  and  hypersonic  external 
flow  show  close  agreement  with  experimental  data  for  both  base  configura¬ 
tions.  Because  of  the  lack  of  clarity  of  the  experimental  exhaust  bound¬ 
aries,  taken  from  photographs,  a  comparison  cotild  only  be  made  to  6  or 
7  nozzle  exit  diameters  downstream  of  the  nozzle  exit. 


Cv:;jLiipeti*ictuxis  of  tlio  supci'aouio  axid  Ixjrpci'ooiiio  nictlxodo  -with,  cxpcri 

mental  boundaries  for  undisturbed  free-stream  Mach  mmibers  of  3.5,  4.0, 
and  5.  0  are  in  good  agreement.  For  cases  with  low  supersonic  and  high 
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hypersonic  external  conditions  a  slight  deviation  of  the  two  methods 
resulted.  Since  no  experimental  boundaries  in  the  low  supersonic  or 
high  hypersonic  free-stream  range  were  available  for  comparisons,  no 
definite  conclusions  could  be  made  as  to  an  upper  or  lower  free-stream 
Mach  number  limit  for  either  method.  However,  based  on  the  results 
presented,  either  method  should  give  a  good  approximation  of  jet  bound¬ 
aries  from  highly  underexpanded  supersonic  nozzles. 

It  is  interesting  to  note  the  effects  of  varying  jet  and  free-stream 
conditions  for  boundaries  calculated  by  the  present  methods.  For  con¬ 
stant  nozzle  geometry,  chamber  total  pressure,  and  altitude,  the  jet 
spreading  increases  with  decreasing  specific  heat  ratio.  Whereas,  for 
constant  specific  heat  ratio,  the  spreading  increases  with  increasing 
altitude  (increasing  Pj/Pa,)  and  decreases  with  increasing  nozzle  exit 
Mach  number  (increasing  Pc/pj)- 
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APPENDIX 

COMPUTATIONAL  PROCEDURE -QUIESCENT  ATMOSPHERE 


Step  1 

a.  Calculate  the  initial  angle  of  jet  boundary  at  the  nozzle  lip  for  a 
base  pressure.  Pb^  =  P„ 


and  vj  correspond  to  Pb^/Pc  ^j»  respectively,  and  are  ob¬ 
tained  from  compressible  flow  tables  based  on  yj.  6]>j  is  known  from 
the  nozzle  geometry. 


b. 


Locate  Point  2  by  assuming  one -dimensional  flow  from  the  nozzle  lip. 
Point  1,  to  some  arbitrary  Point  2  where  from  Eq.  (2) 

Jj-  =  1  +  Ai. 

and  from  Eq.  (3) 


Xj  Ax  A r  cot  Oj 

TT  "  ^  ^ 


step  2 

a.  To  locate  Point  3,  find  the  average  internal  flow  conditions  between 
Points  1  and  2  and  compress  this  flow  to  the  free -stream  static  pres¬ 
sure  at  Point  2.  The  average  internal  flow  conditions  are  calculated 
as  follows: 


1.  Find  (A/A*)i  corresponding  to  p^^/Pc  from  the  compressible  flow 
tables  representing  the  exhaust  gas. 

2.  Calculate  (A/A*)2  from  Eq.  (6)  as 


3.  Find  Mb^  and  Mb2  corresponding  to  {A/A*)^  and  (A/A*)2,  respec¬ 
tively,  from  the  compressible  flow  tables  for  the  y  of  the  exhaust 
gas. 


4.  Calculate  (Mb)avg  between  Points  1  and  2  by  assuming  a  linear 

Mach  number  distribution  over  the  finite  distance  between  1  and  2. 

TP-rnm  (fi), 

Mb.  +  Mbj 


(Mb) 


avg 
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5.  Find  (Pb/Pc^avg  corresponding  to  (Mb)avg  the 
compressible  flow  tables  representing  7j, 

b.  From  Eq.  (4a),  o2  is  obtained  as 


ai  is  known  from  Step  1. 

c.  With  02  determined.  Point  3  may  be  located  by  assuming  one- 

dimensional  flow  from  Points  2  to  3  (same  operation  as  in  Step  lb 
using  Eqs.  (2)  and  (3)). 

Subsequent  points  along  the  boundary  may  be  located  by  repeating 
the  operation  in  Step  2  for  each  point. 


CALCULATIONAL  PROCEDURE  -  HYPERSONIC  EXTERNAL  FLOW 


CONFIGURATION  1 


Step  1 

a.  For  different  values  of  Ofi  calculate  vi  from  the  following  equation: 

=  ai  +  i/j  -  0JJ 

uj  is  obtained  from  compressible  flow  tables  with  Mj  and  yj  for  the 
combustion  gases.  0n  is  known  from  the  nozzle  geometry. 

b.  Find  P^j^/Pc  corresponding  to  »/i  from  the  compressible  flow  tables 
for  the  y  of  the  exhaust  gas. 


c. 

d. 


Plot  Pbi/Pc  versus  oi  from  Steps  a  and  b. 

For  the  same  values  of  ai  calculate  Pbj/Pc  using  Eq.  (7); 


p 

_ 2 

Pe 


fy  M  *  sin*  a  +  1  1 
L  oo  oo  ^  J 


e.  Plot  Pbi/Pc  versus  ai  from  Step  d  on  the  same  curve  as  Step  c. 

f .  From  the  intersection  of  the  two  curves  read  the  values  of  ai  and 
Pbi/Pc- 

g.  Locate  Point  2  by  assuming  one -dimensional  flow  from  the  nozzle 
lip.  Point  1,  to  some  arbitrary  Point  2  where 


Xj  Ax  Ar  cot  Ui 
7^  ^  ?] 
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Step  2 


a.  To  locate  Point  3,  find  the  average  internal  flow  conditions  between 
Points  1  and  2  and  equate  this  flow  with  the  external  flow  as  defined 
by  Eq.  (7).  The  calculations  to  determine  the  internal  flow  condi¬ 
tions  are  as  follows: 


1.  Find  (A/A*)i  corresponding  to  Pb^/Pj,  from  the  compressible 
flow  tables  representing  the  exhaust  gas. 

2.  Calculate  (A/A*)2  from 


3.  Find  Mbi  and  Mb2  corresponding  to  (A/A*)i  and  (A/A*)2.  respec¬ 
tively,  from  the  compressible  flow  tables  for  the  y  of  the  exhaust 
gas. 

4.  Calculate  (Mb)avg  between  Points  1  and  2  by  assuming  a  linear 
Mach  number  distribution  over  the  finite  distance  between  1  and  2 


(Mb) 

avg 


Mbj  +  Mbj 
2 


5.  Find  (Pb/Pc^avg  corresponding  to  (Mb)avg  f^rom  the  compressible 
flow  tables  representing  the  exhaust  gas. 


b. 


For  the  internal  flow,  assume  different  values  of  a2  and  calculate 
Pb2/Pc  using  Eq.  (4): 


(Pb) 


avg 


[yj  (Mb) 


avg 


Sin 


(a,  -  Oj)  +  ij 


c.  Plot  Pb2/Pc  versus  q'2  from  Step  Id  on  the  same  curve  as  defined  by 
the  external  impact  equation  in  Step  Ic. 

d.  From  the  intersection  of  the  two  curves  read  the  values  of  a'2  and 
Pb2/Pc' 

e.  With  Q'2  determined.  Point  3  may  be  located  by  assuming  one¬ 
dimensional  flow  from  Points  2  and  3  (seime  operation  as  in  Step  Ig 
using  Eqs.  (2)  and  (3)). 


Subsequent  points  along  the  boundary  may  be  located  by  repeating 
the  operation  in  Step  2  for  each  point. 
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CONFIGURATION  2 


Step  1 

Calculate  the  constant  pressure  portion  of  the  boundary  (nozzle  lip, 
Point  1,  to  the  point  of  intersection  with  the  free  stream.  Point  i)  as 
outlined  for  the  quiescent  atmosphere. 


Step  2 

a.  Determine  the  internal  flow  conditions  prior  to  interaction  with  the 
free  stream  as  follows:  (Note:  n  =  i  -  1.) 


1.  Find  (A/A*)j  _  j  corresponding  to  p^^  _  j^/Pq  from  the  compres¬ 
sible  flow  tables  representing  the  exhaust  gas. 

2.  Calculate  (A/A*)^  from 


3.  Find  Mbj  _  j  and  Mbj  corresponding  to  (A/A*)^  _  j  and  (A/A*)^, 
respectively,  from  the  compressible  flow  tables  for  the  y  of 
the  exhaust  gas. 


4. 


Calculate  (Mb)g^yg  between  (i  -  1)  and  (i)  by  assuming  a  linear 

Mach  number  distribution  over  the  finite  distance  between  (i  -  1) 
and  (i): 


(Mb) 

avg 


(Mb)._^  +  (Mb). 


5.  Find  (Pb/Pc^avg  corresponding  to  (Mb)avg  from  the  compressible 
flow  tables  representing  yj. 

b.  With  the  average  internal  flow  conditions  between  (i  -  1)  and  (i) 

defined,  the  calculations  for  the  boundary  are  the  same  as  for  Con¬ 
figuration  1,  Steps  2b- 2e. 


CALCULATIONAL  PROCEDURE -SUPERSONIC  EXTERNAL  FLOW 


CONFIGURATION  1 
Step  1 

a.  For  different  values  of  ai  calculate  ui  from  the  following  equation: 

-i  -  -i  '  "j 

uj  is  obtained  from  compressible  flow  tables  with  Mj  and  yj  for  the 
exhaust  gas.  0N  is  known  from  the  nozzle  geometry. 
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b.  Find  Pbi/Pc  corresponding  to  from  the  compressible  flow  tables 
for  the  y  of  the  exhaust  gas. 

c.  Plot  Pb^/Pc  versus  from  Steps  a  and  b. 

d.  For  the  same  values  of  ai,  find  the  corresponding  pressure  coeffi¬ 
cient  Cp  from  conical  flow  charts  based  on  yg^.  The  pressure  coef¬ 
ficient  is  defined  as 

Cp  - 

<ioc  Xoe 

where  Cp  =  f(ai,  M.,). 

e.  Solve  for  Pb^/Pc  follows: 

Pbj  Poo 
Pc  Pc 

y*,  and  p^  are  known  from  the  external  flow  conditions. 

f .  Plot  Pbi/Pc  versus  ai,  from  Steps  d  and  e  on  the  same  curve  as 
Step  c. 

g.  From  the  intersection  of  the  two  curves  read  the  values  of  ori  and 
Pbi/Pc 

h.  Locate  Point  2  by  assuming  one-dimensional  flow  from  the  nozzle 
lip.  Point  1,  to  some  arbitrary  Point  2  where 

=  1  +  Al. 

Xj  Ax  Ar  cot  ai 


step  2 


Successive  points  are  established  by  finding  an  estimated  average 
flow  between  two  preceding  points  and  predicting  the  pressure  at  the  end 
of  the  interval  by  Newtonian  impact  theory.  Equilibrium  is  sought  with 
the  external  flow  which  is  assumed  to  be  governed  by  Prandtl-Meyer 
expansion  of  the  external  flow  behind  the  conical  shock  referred  to  in 
Step  1. 

a.  Define  the  external  flow  aft  of  the  external  shock  as  follows: 

1  .  TTinH  tVio  ov+oT*nQl  TV/TaoVi  niiTr>V>or»,  KoViinrl  +Vio  civ+OT>nal  cKciok 

corresponding  to  and  ai  from  conical  flow  tables. 
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2.  From  compressible  flow  tables  based  on  find  corre¬ 
sponding  to  Mg,  . 

b.  For  different  values  of  Q'2»  calculate  Va,2  follows: 

-  “2) 

c.  Find  Pb2/Pt  corresponding  to  v^2  compressible  flow  tables. 

Pt  is  the  stagnation  pressure  of  the  external  flow  behind  the  conical 
shock. 

d.  Calculate  Pb2/Pc  follows: 

Pc  Pt  V  Pc  /  \  PhJ 

Pbf/Pc  is  obtained  from  Step  Ig  and  p^/Pbi  corresponds  to  v  in 
Step  2b. 

e.  Plot  Pb2/Pc  versus  02  from  Steps  2b  and  2d. 

The  computations  of  the  internal  flow  for  the  step-by-step  equilibrium  of 
internal  and  external  pressures  to  determine  the  jet  boundary  is  outlined 
in  the  following  steps  to  determine  Point  3. 

f .  Calculate  the  average  flow  conditions  between  Points  1  and  2  as 
follows : 

1. 

2. 


3. 


4. 


5. 


Find  (A/A*)2  corresponding  to  Pb^/Pc  from  the  compressible 
flow  tables  representing  the  exhaust  gas. 

Calculate  (A/A*)2  from 

1  '  *  '  2 

Find  Mbi  and  Mb2  corresponding  to  {A/A*)^  and  (A/A*)2,  respec¬ 
tively,  from  the  compressible  flow  tables  for  the  y  of  the  exhaust 
gas. 

Calculate  (Mb)QVfT  between  Points  1  and  2  by  assuming  a  linear 
avg 

Mach  number  distribution  over  the  finite  distance  between  1  and  2. 

(Mb)  =  - - — 


Find  (pb/pc)avg  corresponding  to  (Mb)ayg  from  the  compressible 
flow  tables  representing  the  exhaust  gas. 
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g.  For  the  internal  flow,  assume  different  values  of  a2  and  calculate 
Pb2/Pc  losing  Eq.  (4): 

[yj  (Mb)*  sin"  {a^  -  a^)  +  ll 
Pc  Pc  L  J  avg  J 

h.  Plot  Pb2/Pc  versus  02  from  Step  d  on  the  same  curve  as  defined  by 
the  external  impact  equation  in  Step  2e. 

i .  From  the  intersection  of  the  two  curves  read  the  values  of  a 2 

Pb2^Pc‘ 

3 .  With  02  determined.  Point  3  may  be  located  by  assuming  one¬ 
dimensional  flow  from  2  to  3  (same  operation  as  in  Step  Ih  using 
Eqs.  (2)  and  (3) ). 

Subsequent  points  along  the  boundary  may  be  located  by  repeating  the 
operation  in  Step  2  for  each  point. 

CONFIGURATION  2 
Step  1 

Calculate  the  constant  pressure  portion  of  the  boundary  (nozzle  lip. 
Point  1,  to  the  point  of  intersection  with  the  free  stream.  Point  i)  as  out¬ 
lined  for  the  quiescent  atmosphere. 


Step  2 


a.  Determine  the  internal  flow  conditions  prior  to  interaction  with  the 
free  stream  as  follows:  (Note:  n  =  i  -  1. ) 

1.  Find  (A/A*)j[  _  i  corresponding  to  Pbj  _  ^/Pc  from  the  compres¬ 
sible  flow  tables  representing  the  exhaust  gas. 

2.  Calculate  (A/A*)^  from 


3.  Find  Mbj  _  i  and  Mb^  corresponding  to  {A/A^)^  _  1  and  (A/A*)i, 
respectively,  from  the  compressible  flow  tables  for  the  y  of 
the  exhaust  gas. 


4. 


Calcvdate  (Mb)avg  between  (i  -  1)  and  (i)  by  assuming  a  linear 

Mach  number  distribution  over  the  finite  distance  between  (i  -  1) 
and  (i): 


(Mb) 

avg 


(Mb)._^  +  Mbj 
2 
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5.  Find  (Pfa/Pc^avg  corresponding  to  (Mb)avg  frona  the  compres¬ 
sible  flow  tables  representing  yj. 

b=  Assume  different  values  of  ai  and  calculate  Pbi/Pc  using  Eq.  (4): 


Pbi  _  (Pb) 


avg 


[yj  (Mb) 


Sin 


avg 


(ai_i  -  Oi)  +  ij 


c.  Plot  Pbi/Pc  versus  ai  from  Step  2b. 

d.  For  the  same  values  of  ori,  find  Pbi/Pc  the  external  flow  from 
Steps  Id  and  le  for  Configuration  1. 

e.  Plot  Pbi/Pc  versus  ai  from  Step  2c  on  the  same  curve  as  Step  2d. 

f .  From  the  intersection  of  the  two  curves  read  the  values  of  ai  and 

Pbi/Pc- 

g.  Locate  Point  (i+  1)  as  in  Step  Ih  for  Configuration  1. 

Successive  points  are  obtained  by  repeating  the  operations  outlined 
in  Step  2  for  Configuration  1. 
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Jet  Boundary 

Still  Air 


External  Flow 


Jet  Boundary 
External  Flow 
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Boundary  Static  to  Chamber  Stagnation  Pressiure  Ratio 
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Fig.  3  Initial  AngU  Plot 
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O.  D/D|  .  1.0,  M|  -  2.5 


b.  D/D;  o  2.0,  M,  -  3.7 


Note:  All  dimensions 
in  inches 


Fig.  4  Schematic  of  Toot  Modeli 
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CD 


Boundary  Mixing  Region 


Present  Boundary  Solution 


Jet  Shock 


Fig,  5  Typical  Schlieren  Photograph  of  the  Jet  Structure,  =  5.0 
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to  Exit  Radius  Ratio,  r/r 


Fig*  6  Jet  Boundaries  Calculated  by  the  Present  Quiescent  Method  Compared  with  the  Method  of  Characteristics  Solution 
and  Latvala's  Solution  for  a  Quiescent  Atmosphere 
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Fig.  7  Jot  Boondarloo  Calculatod  by  ifao  ProEont  Hyportontc  Molhod  Comparod  with  tho  Method  of  Characteristics  Solntion 
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Axial  Distance  to  Ex 
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Radial  DisLaxice  to  Exit  Radius  Ratio,  r/r 


Fig.  8  Jet  Boundaries  Calculated  by  the  Present  Hypersonic  Method 
Compared  with  Experimental  Jet  Shocks 
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b.  =  5.0 


Fig.  8  Concluded 
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Radial  Distanr-u  to  Exit  Radius  Ratio,  r/ 


0  2  4  6  8  10  12  14 


Axial  Distance  to  Exit  Radius  Ratio,  x/rj 

a.  Pj 10,  Pc/Pj  =  101 


Fig.  9  Jet  Boundaries  Calculated  by  the  Present  Methods  Compared  with 
Experimental  Boundaries,  -  3.5 
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Radial  Distance  to  Exit  Radius  Ratio,  r/r^ 


1.4  D/Dj  =  2 


Experimental  Mixing  Boundar^y 


fj  =  17  deg 


Calculated  Jet  Shock 


- Experimental  Jet  Shock 

O  Hypersonic  Boundary  Solution 
n  Supersonic  Boundary  Solution 


Axial  Distance  to  Exit  Radius  Ratio,  x/rj 
b.  pj/p_^=  15,  p^/pj=  101 

Fig.  9  Concluded 
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Radial  Distance  to  Exit  Radius  Ratio,  r/r ^ 
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Axial  Distance  to  Exit  Radius  Ratio,  x/r  . 
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Fig.  10  Jet  Boundaries  Calculated  by  the  Present  Methods  Compared  with 
Experimental  Boundaries,  M  >  4.0 
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Radial  Distance  to  Exit  Radius  Ratio,  r/r-i 
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Experimental  Mixing  Boundary 
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Axial  Distance  to  Exit  Radius  Ratio,  x/r ^ 

b.  p./p^  =  20.0,  p,/pj  =  101 


Fig.  10  Concluded 
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Radial  Distance  to  Exit  Radius  Ratio,  r/r 


Axial  Distance  to  Exit  Radius  Ratio,  x/rj 

«•  Pj/P„  =  35,  Pc/pj  =  101 


Fig<  11  Jet  Boundaries  Calculated  by  the  Present  Methods  Compared  with 
Experimental  Boundaries,  M^  ^  5,0 
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Radial  Distance  to  Exit  Radius  Ratio,  r/r^ 
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Radial  Distance  to  Exit  Radius  Ratio,  r/rj 


Axial  Distance  to  Exit  Radius  Ratio,  x/rj 


a.  Configuration  1 
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Fig.  13  Coitelndad 
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Radial  Distance  to  Exit  Radius  Ratio,  r/rj 


Fig.  14  Jet  Beendaries  Calculated  by  the  Present  Hypetsenic  Method  lor  Various  Nezsie  Mach  Nusibers 
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